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Ah initio calculations have been carried out to study the magnetic dipole and electric 

^ , quadrupole hyperfine structure constants of ^*'^Pb+. Many-body effects have been con- 

^ ■ 

CJ I sidered to all orders using the relativistic coupled-cluster theory in the singles, doubles and 

Z^ ' partial triples approximation. The trends of these effects are found to be different from 

1-5 ' atomic systems that have been studied earlier. 

>^ I. INTRODUCTION 

oo . 

With the advent of ion trapping techniques, it has become possible to perform high precision 
<^ ' measurements of different atomic properties; including hyperfine structure constants for the ground 



^: ' ' riQ 

O . and low-lying excited states of stable isotopes |lj, |2| 



c/3 
Ph. 



_ ^ _^ Studies of these interactions have served as 
Q ' stringent tests of relativistic many-body theories. It has been found that the behavior of electron 



correlation in the hyperfine interactions in the d-states is substantially different from those of s- 
and p- states |3]. The underlying reason for this is the strong and unusual core-polarization effects 

associated with the former. 

• 1—1 , 

r> \ The work presented in this paper is carried out using the relativistic coupled-cluster (RCC) 

theory. This theory has certain features which make it suitable for the calculations of excited 
state properties. Prominent among them being its abilities to treat relativistic and correlation 

n 

effects in the initial and final states in a rigorous and balanced way '4\. It is equivalent to all 
orders relativistic many-body perturbation theory. Core-polarization and pair-correlation effects 
which are important for such calculations are computed to all orders in this theory. It has been 
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successfully applied to heavy atomic systems with a single valence electron |J, |5|| . 

205pj-, jg Qj-^g q£ ^YiQ stable isotopes of lead which has lifetime about 1.5 x 10^ years. The nuclear 
structure of this isotope is different from the other stable isotope, ^'^^Pb. It's nuclear magnetic 
moment and electric quadrupole moment are non-zero. Theoretical studies of the magnetic dipole 
and electric quadrupole hyperfine structure constants of 205p^+ ^^.g ^^le focus of the present work. 

There are relatively few measurements of the hyperfine structure constants for the excited states 
of heavy atoms and ions. It is certainly worthwhile to perform highly correlated calculations of 
these quantities for Pb'^ as it could be a promising candidate for the observation of parity non- 
conservation ja]. This theoretical effort could motivate experimentalists to carry out high precision 
measurements of ^'^^Pb^. 

In section II of this paper, we give an outline of the RCC theory and in section III we discuss 
the pertinent details of the calculation. The results of our calculations are presented and discussed 
in section III and finally we make some concluding remarks in section IV. 

II. THEORY 

In order to obtain the RCC wavefunction for Pb"^, we require the closed-shell Pb++ wavefunc- 
tion, which is given in terms of the lowest order Dirac-Fock (DF) wavefunction \^£)f), as 

\^) = n\^DF), (2.1) 

where in conventional many-body perturbation theory, the wave operator, 0, is expressed in powers 
of the residual interaction, Ves = J2i<j ^t' ~ ^i ^DF{i"i)- This results in a hierarchy of approxima- 
tions for the correlation energy and the wavefunction. 

The coupled cluster theory is based on the following exponential ansatz for the closed-shell wave 
operator [TJ 

n = N[exp{T)], (2.2) 

where the cluster operator T is expressed in terms of the connected diagrams of the wave operator. 
The operator T also can be expressed in terms of the cluster operators T„ corresponding to different 
orders of excitations n of the core electrons from the DF state, \^df), explicitly defined by, 

T = T1 + T2 + ... 

ap abpq 



with a,b,c, ..{p,q,r, ..) representing occupied (unoccupied) orbitals. t^ is the cluster ampHtude 
corresponding to the single excitation from the orbital a to p and so on. Termination of the 
series at T2, results in the coupled cluster J;heory with single and double excitations (CCSD). The 
contribution from the Breit interaction 8^ which is four orders of magnitude smaller than the 
Coulomb interaction has been neglected in the present work. 

For a single valence atomic system the wavefunction in the RCC method can be written as 



with the new reference state 



|^„) = e^{l + 54|$„), (2.4) 



|$.)=at|$BF), (2.5) 



for the given valence electron v; S^ represents excitation operators which excite at least the valence 
electron. The explicit form of this operator can be written as 

Sv = Sly + S2V + ••• 

Py^v bpq 

An approximate treatment of the triple excitations to the CCSD method is included by con- 
tracting the residual Coulomb operator, which effectively forms a two-body operator, and the 
double-excitation operators T2 and S2V lid], thereby defining the CCSD(T) approximation 

tij ~r fcfo ~r tc tp tq t^ 

where e^ is the orbital energy of the corresponding i 'th electron. 



III. METHOD OF CALCULATION 

For computational simplicity, the T amplitudes are solved first for the closed-shell Pb++ and 
then the valence electron is attached to calculate the open-shell wavefunctions. The matrix equa- 
tions for solving the correlation energy and the T amplitudes are given by 

{^DfI'^I^Df) = AEcorr (3.1) 

i'^hpi^l'^DF) = 0, (3.2) 

where Hjv is the normal ordering Hamiltonian which can be written as Hjv = f at + Vat, sum of 
one-body and two-body terms. In the above equation. Hat is defined as Hat = e~-^HAre^, AEcorr 
is the correlation energy and \^*£)p) corresponds to excited states from \^£)p). 



For the open-shell RCC calculations, the S^j operators are solved using the following equations 

($„|H^{1 + 54|$„) = -AE, (3.3) 

($:|H^{1 + S4I$,) = -AE,{^:\{S,}\^,), (3.4) 

where AE^ is the ionization potential (IP) energy of the corresponding valence electron v. Similarly 
!<!>*) are the excited states from \^v)- 

The relativistic hyperfine Hamiltonian is given by 11 1 



if;,/, = ^M«.tW, (3.5) 

k 

where M''^^ and T''^^ are spherical tensor operators of rank k. In first-order perturbation theory, 
hyperfine energies E^fsiJ) of the fine-structure state \JMj) are expectation values of the hyperfine 



interaction Hamiltonian. Details of the expression are given by Chang and Childs 
magnetic dipole and electric quadrupole hyperfine energies are defined by 






The 



Emi = AK/2 (3.6) 

and 

_g 3i^(i^ + l)-4/(/ + l)J(J + l) 
^^ 2 2/(2/- 1)2 J(2J-1) ' ^ ' 

respectively. Here / and J being the total angular momentum of the nucleus and the electron state, 

respectively, and K = 2{I ■ J). The magnetic dipole hyperfine constant A and electric quadrupole 

hyperfine constant B are defined as 

(J||tW||J) 

A = LlNgi^l= , (-J-o) 

^ " ^J{J+l){J+l) 

and 

2J(2J-1) T^/^ 



B = 2eQ 



(j||r(2)||j), (3.9) 



X2J + l){2J + 2){2J + 3)_ 
respectively, where fi^- is Bohr magneton, gj = fii/I with m and / are the nuclear dipole moment 
and spin, and Q is nuclear quadrupole moment. 

The expectation value for a general one particle operator O in a given valence electron (v) state 
can be expressed in RCC theory as 



(O). 



{<^,\{l + Sl}0{l + S,}\<P,) 



where we define 



and 



= e^^Oe^ (3.11) 



= ($,|5K+^K^. + ^'5il^^-)- (3.12) 

For computational simplicity we evaluate the matrix elements of any operator in two steps. We 
expand O using Wick's general theorem |7| as 

O = (e^'Oe^)/.e. + (e^'Oe^)„.b. + {e'^' Oe^U. + ...., (3.13) 

where we have used the abbreviations f.c, o.b. and t.b. for fully contracted, effective one-body 
and effective two-body terms respectively. In this expansion of O, the effective one-body terms are 
computed keeping terms of the form of 

Oo.b. = O + T^O + OT + T^OT. (3.14) 

The calculation procedure for these terms are given by Geetha et al ja|. They are finally 
connected with S^ and Sl operators in the evaluation of properties . Contributions due to the 
effective two-body terms from O are constructed using the procedure shown diagrammatically in 
our earlier works [a, llj] and computed directly during the calculation of properties. The following 
The following types of terms are considered for the construction of the effective two-body terms 

Ot.b. = OTi + tIo + OT2 + T^O. (3.15) 

Other effective terms correspond to higher orders in the residual Coulomb interaction and are 
neglected in the present calculation. A similar procedure has been followed to account for the 
normalization factor. 

The contributions from the normalization factors for the corresponding valence electron v that 
are given in in tables III and IV, were obtained using the following relations: 

Norm={^^\0\^y){-—---l}. (3.16) 

J- + ^^v 

IV. RESULTS AND DISCUSSIONS 

The starting point of our calculations is the generation of DF orbitals for ^'^^Pb^^. These 
orbitals are constructed as linear combinations of Gaussian type orbitals (GTOs) as given by 



Chaudhuri et al. [ij]. It has been found |a| that the RCC calculations based on these kind of 
GTOs provide accurate A— values for the low- lying states of ^'^ Pb+. In the present work, we 
have used similar basis functions to calculate both A— and B— values of ground and some of the 
important excited states of ^'^^Pb"'". 

In table I, we present A and B hyperfine structure constants of the low-lying states of Pb"^. We 
use the Lande nuclear g- fa ct or, gj = 0.28468 to calculate A and nuclear quadrupole moment, Q 
= 0.234 to calculate B [lj|. From the differences of the DF and RCC results given in table I, it 
is evident that for both the hyperfine constants A and B the electron correlation effects vary from 
(10-290)%. 

TABLE I: A and B results of 205pb+ using DF and RCC methods. 





6pi/2 


6P3/2 


7si/2 


7pi/2 


7P3/2 


Qd-i/2 


6^5/2 


8si/2 


A 


















DF 


2765.54 


220.58 


1879.06 


476.34 


43.84 


19.05 


7.53 


687.77 


RCC 


3099.5 


149.7 


2680.3 


543.7 


74.8 


-10.03 


60.97 


929.3 


B 


















DF 




377.81 






75.09 


11.08 


13.07 




RCC 




464.6 






99.3 


50.7 


56.9 





All the core orbitals were excited in our calculations. The core-polarization effects, which are 
the largest contributors to the hyperfine constants of the 6^3/2 state of ^'^^Pb"'" 6] and the d^/2 
states of the alkaline earth ions [ISl, have been accounted to all orders through the 05*2?; term. 
It was also found from the hyperfine structure studies of the Si /2 and pi/2 states in the alkaline 
earth ions that pair-correlation and core-polarization effects are important. In order to appreciate 
the importance of these effects in Pb^, we present their contributions in table II. Comparison of 
these results with their corresponding DF values from Table I, brings out some distinct many-body 
features of the system. The most prominent among them is the size of the core-correlation and 
core-polarization effects for the d— states. As in the case of some of the alkaline earth ions, the 
sign of the core-polarization effect in the ^5/2 state is opposite that of the DF value and the net 
contribution is 209% of the corresponding DF value. However, unlike those alkaline earth ions 
[l5| the final RCC result has the same sign as the DF result. This is due to the large positive 



TABLE II: The contributions of core-correlation^ core-polarization and pair-correlation of A and B results 
in 



205p]-,+ 



Statos Core-corr. Core-pol. Pair-corr. 

A 

7p3/2 -0.68 15.17 6.01 

6^3/2 -1.02 11.13 9.30 

6^5/2 -0.48 -15.73 3.45 

8s 70.29 144.23 139.92 

B 

6p3/2 29.26 79.19 35.55 

7p3/2 3.34 15.78 10.29 

64/2 -0.65 17.87 5.44 

6^5/2 -0.88 23.46 5.96 



contributions from the S2yOS2v term, shown in table III . Therefore, it is imperative to consider 
such terms to obtain accurate results. 

The role of electron correlation in the hyperfine interactions in '^^Sc and ^^Y has similarities 
with Pb"*" even though those systems are neutral and have different electronic configurations. 
Another interesting feature of the present study on 205p)-,-i- jg ^j^^t the RCC result of the A— value 
of 6^3/2 is of opposite in sign than that of the DF result. This trend is different from the other ^3/2 
states in alkaline earth ions l3, llol- The main reason for this behavior is due to another unusual 



contribution from S2yOS2v The core correlation effect on B— values though reduces as higher 
excited states are considered, which is expected, but contributions to the d— states are not similar 
to the p— states. This is evident from the B— values given in table II. 

Table III and IV present the important effective two-body terms obtained from O = e^ Oe^ , 
but they contribute very little. One can therefore justifiably ignore the higher order terms given in 
eqn. (3.15) and save computational time. The correlation effects of the A— values of other states 
presented in table I, behave the same way as in ^'^'^Pb'*' 



TABLE III: Contributions of different coupled-cluster terms to the Pb+ magnetic dipole hyperfine structure 
constant (A), cc stands for the complex conjugate part of the corresponding terms. 



Terms 


7p3/2 


6d3/2 


6d5/2 


8Si/2 




state 


state 


state 


state 


(DF) 


43.84 


19.05 


7.53 


687.77 





44.52 


20.07 


8.01 


617.48 


OSiy + cc 


6.01 


9.30 


3.45 


139.92 


OS2V + cc 


15.17 


11.13 


-15.73 


144.23 


SlOSi, 


0.21 


1.09 


0.38 


7.92 


Sl^0S2v + cc -0.86 


2.02 


-2.08 


4.85 


SlJ)S2v + cc 11.04 


-54.93 


68.91 


29.74 


Important 


effective two-body terms of 



Sl^OTi+cc -0.03 -0.03 -0.01 1.03 

S\^OT2 + cc -0.44 0.74 0.26 -3.61 

Norm -0.69 0.39 -2.21 -10.06 

V. CONCLUSION 

The RCC theory has been employed to study the magnetic dipole and electric quadrupole 
hyperfine structure constants of the ^'^^Pb^. Strong electron correlations effects are found in the 
d— states and their behavior is different from other systems studied earlier. Experiments to measure 
these quantities will constitute important tests of the relativistic coupled-cluster theory. 

VI. ACKNOWLEDGMENT 

We are grateful to Prof. Werth for valuable discussions and suggestions for this calculation. 
The calculations were carried out using the Tera-flop Supercomputer at C-DAC, Bangalore, India. 



[1] X. Feng, G.-Z. Li, R. Alheit and G. Worth, Phys. Rev. A, 46, 327 (1992) 

[2] W. Liu et al, Phys. Rev. Lett, 82, 711 (1999) 

[3] B.K. Sahoo et al., Phys. Rev. A, 68, 040501 (R) (2003); B.K. Sahoo et al, J. Phys. B,38, 4379 (2005) 



TABLE IV: Contributions of different coupled-cluster terms to the Pb+ electric quadrupolc hypcrfine struc- 
ture constant (B). 



Terms 6p3/2 7p3/2 edg/s Odg/a 

state state state state 

O (DF) 377.81 75.09 11.08 13.07 

O 348.55 71.56 11.73 13.95 

OSiy+cc 35.55 10.29 5.44 5.96 

0S2V+CC 79.19 15.78 17.87 23.46 

SlJJSi^ 0.87 0.36 0.64 0.65 

SlJ)S2v + cc 2.18 0.62 2.27 2.75 

sIjJS2v + cc 10.11 2.56 14.26 11.75 

Important effective two-body terms of O 

sLOTi + cc -0.72 -0.05 -0.02 -0.02 
Sl,0T2 + cc -5.19 -0.76 0.43 0.46 

Norm -4.89 -0.93 -2.02 -2.07 

[4] B. K. Sahoo, S. Majumder, R. K. Chaudhuri, B. P. Das and D. Mukherjee, J. Phys. B, 37, 3409 (2004) 
[5] Geetha Gopakumar, Holger Merlitz, Rajat Chaudhuri, B.P. Das, U.S. Mahapatra, and D. Mukherjee 

66, 032505 (2002) 
[6] B. K. Sahoo, R. K. Chaudhuri, B.P. Das, H. Merlitz and D. Mukherjee, Phys. Rev. A 72, 032507 (2005) 
[7] I. Lindgren and J. Morrison, Atomic Many-body Theory 3, ed. G. E. Lambropoulos and H. Walther 

(Berlin: Springer) (1985) 
[8] G. Breit, Phys. Rev. 34, 553 (1929); 34, 383 (1930); 39, 616 (1932) 
[9] D. Mukherjee and S. Pal, Adv. Quantum Chem. 20, 281 (1989) 032507 (2005). 
[10] A. Haque and U Kaldor, Chem. Phys. Lett. 117, 374; 120, 261 (1985) 
[11] K.T. Cheng and W.J. Childs, Phys. Rev. A, 31, 2775 (1985) 
[12] B. K. Sahoo, PhD thesis submitted to Mangalore University, India (2005) 
[13] R. K. Chaudhuri, P. K. Panda, B. P. Das, U. S. Mahapatra, and D. Mukherjee, J. Phys. B 33, 5129 

(2000) 
[14] M. Anselment, W. Faubel, S. Groing, A. Hanser, G. Meisel, H. Rebel and G. Schatzz, Nucl. Phys. A, 

451, 471 (1986) 
[15] B. K. Sahoo, C. Sur, T. Bcicr, B. P. Das, R. K. Choudhuri and D. Mukherjee, submitted to PRL (2006) 
[16] C. Sur, B. K. Sahoo, R. K. Chaudhuri, B. P. Das and D. Mukherjee, Eur. Phys. J. D 32, 25 (2005) 



